Electronic states and optical properties of GaAs/AlAs and GaAs/vacuum 
superlattices by the linear combination of bulk bands method 
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The linear combination of bulk bands method recently introduced by Wang, Franceschetti and 
Zunger [Phys. Rev. Lett. 78, 2819 (1997)] is applied to a calculation of energy bands and optical 
constants of (GaAs)„/(AlAs)„ and (GaAs)„/(vacuum)„ (001) superlattices with n ranging from 4 
to 20. Empirical pseudopotentials are used for the calculation of the bulk energy bands. Quantum- 
confined induced shifts of critical point energies are calculated and are found to be larger for the 
GaAs/vacuum system. The Ei peak in the absorption spectra has a blue shift and splits into 
two peaks for decreasing superlattice period; the E2 transition instead is found to be split for 
large-period GaAs/AlAs superlattices. The band contribution to linear birefringence of GaAs/AlAs 
superlattices is calculated and compared with recent experimental results of Sirenko et al. [Phys. 
Rev. B 60, 8253 (1999)]. The frequency-dependent part reproduces the observed increase with 
decreasing superlattice period, while the calculated zero-frequency birefringence does not account 
for the experimental results and points to the importance of local-field effects. 



PACS numbers; 73.20.Dx, 78.66.Fd, 71.15.Hx, 78.20.Bh 



I. INTRODUCTION 

Quantum confinement in semiconductor heterostruc- 
tures modifies the energy and dimensionality of electronic 
levels and leads to blue shifts of the optical gaps. Most 
experimental investigations have focused on the energy 
region of the fundamental gap, which is easily accessible 
by photoluminescence and photoluminescence excitation 
spectroscopies and yields a variety of interesting physical 
phenomena related to bound excitonic states; relatively 
few studies of confinement effects on high-energy transi- 
tions have been presented. Concerning theory, confined 
electronic levels close to band edges and the resulting op- 
tical properties can be calculated rather simply and accu- 
rately by the envelope-function method. The theoretical 
problem of determining optical spectra of semiconductor 
heterostructures in the whole visible region is much more 
complex and beyond the reach of effective-mass methods, 
as it requires a description of the effects of confinement 
on electronic states in the whole Brillouin zone. 

Interband absorption spectra of tetrahedral semicon- 
ductors are domin|ated by two prominent features, de- 
noted El and £'2.EJi3 The Ei peak (and its spin-orbit 
counterpart Ei -I- Ai) originates from band-to-band tran- 
sitions along the T — L direction, where valence and con- 
duction bands are nearly parallel: this results in a Mi- 
type critical point, i.e., a saddle point in the joint density 
of states, which also gives a strong excitonic character to 
the transition. The E2 peak, instead, has contributions 
from different parts of the Brillouin zone, but mainly 
from a region centered around the special point (|, |, |) 
(in units of 27r/a, where a is the lattice constant). The E2 
peak has essentially no excitonic character. Blue shifts 
and splittings of the Ei and E2 transitions were measured 



in GaAs/AlAs superlatticesBa. More recently, a quan- 
tum confinement induced shift of Ei and E2 was nwSr 
sured in Ge nanoparticles embedded in a glassy matrixH^El. 

A main purpose of this work is to calculate the be- 
havior of El and E2 transitions upon confinement in 
GaAs/AlAs superlattices and in free-standing GaAs lay- 
ers, which are simulated by GaAs/vacuum superlattices. 
Besides the obvious interest of the GaAs/ Al As system, 
the motivation for choosing GaAs/vacuum superlattices 
is to study confinement effects on Ei and E2 in a sys- 
tem where the electronic states are truly confined in the 
GaAs layers even at high energies. In GaAs/AlAs super- 
lattices, on the other hand, the band structures of the 
two constituents far away from the fundamental band 
edges are rather similar and strong banding effects oc- 
cur in short-period structures, i.e., the electronic states 
become delocalized along the superlattice. A compari- 
son between the two systems should therefore elucidate 
the respective roles of quantum confinement and super- 
lattice band formation in determining the optical proper- 
ties. Free-standing GaAs films can be produced by chem- 
ical etchingZl; GaAs/vacuum superlattices can also be a 
model for multilayers formed by GaAs and a wide-gap 
oxide, like AI2O3 or oxidized AlAs (AlOx). 

The electronic structure of short-period (GaAs)„/ 
(AlAs)„ superlattices can be calculated from first prin- 
ciples by densitypfunctional theory and norm-conserving 
pseudopotentialsu. The computing, time grows rapidly 
with the number n of monolayers.cl Moreover, density- 
functional calculations j-ef the optical properties suffer 
from the gap problenJ13_.and require the inclusion of 
quasiparticle correctionsEj and excitonic effectsjl3 thus 
increasing the computational cost. Receritiy a new pseu- 
dopotential approach has been developedEj, in which the 
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electronic states of a heterostructure are represented as 
linear combinations of bulk bands (LCBB). This method 
can be applied to heterostructures with a large number 
of atoms per primitive cell with a reasonably low com- 
putational effort; it is particularly suited to study how 
the optical spectra of the bulk are modified upon con- 
finement, since the electronic states of the bulk are the 
starting point of the technique. The LCBB method is 
adopted in this work. The results presented here can be 
the starting point for tackling more complex problems 
(e.g. the inclusion of local-field and cxcitonic effects, or 
the study of confinement in large quantum dots) as well 
as for comparison with ab-initio calculations when they 
will become available. 

While GaAs and AlAs are cubic materials and have 
an isotropic dielectric constant, (OOl)-grown superlattices 
are intrinsically uniaxial and the component e\\ of the 
dielectric tensor parallel to the growth direction differs 
from the perpendicular component ej_- The linear bire- 
fringence An = (e^)^/^ p_,[£i|)^/^ has been measured in 
GaAs/AlAs superlatticesllall3: An is of the order of a few 
percent and has a nontrivial dependence on the superlat- 
tice period. Birefringence is much larger in GaAs/AlOx 
multilayers, where it has been employed tOp^hieve phase 
matching for second-harmonic generationtZI. A soal of 
the present work is to calculate linear birefringencal3 and 
its dependence on energy and superlattice period. The 
calculations presented here yield only the contribution to 
birefringence arising from quantum-confinement-induced 
modifications of the electronic states; they do not ac- 
count for the intrinsic dielectric anisotropy of a multi- 
layer arising from different boundary conditions for ..an 
electric field parallel or perpendicular to the layerstBa. 
This second contribution to birefringence is in fact equiv- 
alent to thjeriHclusion of local-field effects in the dielectric 
responseE^c3 and is left for future work. 

The rest of this paper is organized as follows. In 
Sec. II we briefly describe a few aspects of the adopted 
method: the LCBB technique combined with empirical 
pseudopotentials|13u as well as space symmetry and the 
tetragonal Brillouin zone of the superlattice. In Sec. Ill 
we present the results for GaAs/ AlAs and GaAs/ vacuum 
superlattices. These are divided into three parts: elec- 
tronic levels (with the added complication of surface 
states in the GaAs/ vacuum system), optical properties 
via the real and imaginary parts of the dielectric func- 
tion, and birefringence. Section IV contains a discussion 
of the present findings also in view of extensions of this 
work. 



II. METHOD OF CALCULATION 

We consider (GaAs)„/ (AlAs)„ and (GaAs)„/(vacuum)i 
superlattices grown along a (001) crystallographic sur- 
face and with period d = na with n ranging form 4 to 20. 
The ideal structure for a lattice-matched system with 



abrupt interfaces is a simple tetragonal Bravais lattice, 
with a supercell defined by the basis vectors (1, 1, 0) a/2, 
(— 1, 1, 0) a/2, (0, 0, 1) na, where a is the bulk lattice con- 
stant. The reciprocal lattice is also simple tetragonal, 
with basis vectors (1, 1, 0) (27r)/a, (— 1, 1, 0) (27r)/a and 
(0, 0, 1) (27r)/(na). The first Brillouin zone is shown in 
Fig. 1. Superlattice high symmetry points are distin- 
guished from their bulk counterparts by putting a bar 
over the symbol. An additional symmetry point L is 
defined as follows: L = X if n is even, L = .R if n is odd. 
The most important zinc-blende k points are folded onto 
superlattice points as follows: 

L J j=-n+l,n 

Lni,\-X'\ -^L. (1) 

In the case of a common anion structure like 
(GaAs)„/(AlAs)„ the point group is D2d, otherwise it 
is C21,: the iatter is the case of (GaAs)„/(vacuum)„ 
superlattices.ES 

In ordejj4p study large scale systems, we adopt a LCBB 
methodll^llJ: superlattice electronic wave functions are 
expressed as linear combinations over band indices n^ 
and wave vectors k of full-zone Bloch eigenstates of the 
constituent bulk materials. Because of the requirement 
of periodicity, the superlattice potential mixes only bulk 
states labelled by k vectors which differ by a superlat- 
tice reciprocal lattice vector Gsl- the number of cou- 
pled states is hence always equal to 2n, because exactly 
2n Gsl are contained in the fee Brillouin zone. The 
maximum dimension of the basis set is then given by 2n 
multiplied by the number of selected bulk bands indices. 
In the case of a GaAs/vacuum superlattice we decide to 
include the 4 valence bands and the 4 lowest conduction 
bands in the basis set. In the case of a GaAs/ AlAs super- 
lattice the roughest selection is to take both GaAs and 
AlAs bulk states at each mixed k and n;,, orthonormal- 
izing at the end the basis set obtained. But, as GaAs 
and AlAs band structures are very similar except for the 
lowest conduction band (see Fig. || and details of calcula- 
tion below), we have verified that it is enough to include 
only GaAs states for n^ from 1 to 8 together with the 
5th band of AlAs (i.e., the lowest conduction band). The 
resulting set must be orthonormalized. It can easily be 
seen that the final dimension of the basis is always small 
(40 X 9 for the largest supercell). A convergence test is 
presented at the beginning of Sec. III. 

Unlike tight-binding or standard plane wave expan- 
sions, the LCBB method allows the intuitive pre-selection 
of physically important states to be included in the ba- 
sis set. By contrast with the k • p approach, off-F states 
Un.k^o are directly considered, permitting a correct treat- 
ment of T-X coupling and of all confinement-induced cou- 
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plings within the Brillouin zone, without the need for a 
large basis of k = bulk states. The formalism applied 
here is strain-free, believing that the small lattice mis- 
match (~ 0.15%) justifies a description of the geometrical 
structure as unrelaxed, with a lattice constant of 5.655 
A averaged over the constituent materials. All trends in 
the superlattice states obtained by LCBB method were 
shown to be reproducedllj, with a surprising accuracy 
(10-20 meV) and a small computational effort, down to 
the monolayer superlattice and up to large periods. 

A local empirical pseudopotcntial approach is chosen 
to perform the band structure calculations, both for the 
heterostructures and the constituent compounds. Since 
the adopted pseudopjOtentials are designed for a kinetic- 
energy cutoff of 5 Ry,c3 bulk eigenfunctions are expanded 
on a basis set made of about 60 plane-waves at each k 
point. In Fig. ^ we show the band structures of GaAs 
and AlAs calculated with these pseudopotentials. 

The potential term in the one-particle Hamiltonian 
is built as a superposition of screened, spherical atomic 
pseudopotentials v^- 

E «a(r-R-d„) W^„(R) , (2) 

a Re_DL 

where R is a fee direct lattice (DL) vector and do, the 
displacement of the atom of type a in the bulk primi- 
tive cell. The weight function Wa (R) selects the atom 
basis which lies on each lattice site, describing the geo- 
metrical details and the symmetry of the structure: in 
the vacuum layers its value is zero. The Hamiltonian 
matrix elements on the bulk basis set depend on the 
Fourier transform of the pseudopotentials Va (r); details 
on the LCBB method can be found in Refs. [13|jl4| . We 
use the continuous-space functions Va (k) proposed by 
Mader and Zungeil23. In Ref. |2^ the empirical parame- 
ters are adjusted in order to fit both the measured elec- 
tronic properties of bulk GaAs and AlAs and some local 
density approximation (LDA) results for superlattices. 
It has been verified that the wave functions of bulk and 
superlattices systems calculated with these pseudopoten- 
tials are close to ttic^e obtained in rigorous first principles 
LDA calculationscEl. The As potential depends on the 
type of nearest neighbors (Ga or Al), thereby describing 
the local-environment dependence of the atomic poten- 
tial, by considering the local electronic charge. These 
pseudopotentials are adjusted to reproduce the experi- 
mental GaAs/AlAs valence-band offset (0.50 eV). Bulk 
and superlattice energy levels are provided in the same 
absolute energy scale, thus superlattice and bulk eigen- 
values can be compared directly. The inclusion of non- 
local terms is needed to achieve a better description of 
high energy states, nevertheless we have verified that al- 
most only transitions to the levels originated by the two 
lowest conduction bands are responsible for the struc- 
tures in the optical spectra below 6 eV, therefore non- 
local terms can be neglected for our purposes. At this 
stage we have also decided to neglect spin-orbit interac- 



tion, even if in Ref. |2^ it is suggested how to include it. To 
preserve a correct description of interfaces in GaAs/ AlAs 
superlattices, an As atom bound to two Al and two Ga 
atoms has been attributed a symmetrized pseudopotcn- 
tial which is the average of the As pseudopotcntial func- 
tions in GaAs and AlAs environments. In GaAs/ vacuum 
superlattices, the dangling bonds at the interfaces are re- 
sponsible for the appearance of surface states, lying in the 
forbidden energy gap. The sizeable metallic contribution 
to low energy transitions is an artifact of the calculation 
(in a real system dangling bonds would be saturated by 
the microscopic nature of the interface) and must be min- 
imized by excluding surface states from initial and final 
states in optical transitions. 

Starting from the one-electron band structure, the 
complex dielectric function e{Lu) = ei -I- ie2 is evalu- 
ated in a straightforward way hyjieans of semiclassical 
theory of interband transitions.oO Tetragonal Brillouin 
zone integrations are performed using Fourier quadra- 
ture with 1056 CbFi|di. and Cohen special points in the 
irreducible wedge.Eao An empirical gaussian broaden- 
ing of 0.1 eV has been introduced to obtain smoothed 
curves. Bulk spectra have been preliminarly calculated 
integrating over both a fee Brillouin zone and the corre- 
sponding folded tetragonal Brillouin zone, verifying the 
exact equivalence of the two results. 

In Fig. ^ we show the calculated real and imaginary 
parts of dielectric functions for bulk GaAs and AlAs crys- 
tals. Agreement with the experimental positions of the 
peaksH'cl is within a few tenths of an eV: we underline 
that the empirical approach is free of the band gap prob- 
lem. The height of the peaks, especially Ei, cannot be 
correctly estimated without the inclusion of the excitonic 
contributions and the local field effects. 



III. RESULTS 

A. Electronic levels 

After calculating the single-particle eigenstates of both 
bulk systems in the whole Brillouin zone, in order to ex- 
tract the bulk states basis on which to construct the su- 
perlattice Hamiltonian, we applied the method described 
in section II to obtain the single-particle band structures 
of (001) (GaAs)„/(AlAs)„ and (GaAs)„/(vacuum)„ su- 
perlattices. The period d = na has been varied from 4a 
to 20a. 

When a sufficiently large number of bulk states is used 
as a basis set for the LCBB method, the results must con- 
verge to those obtained with a direct diagonalization of 
the Hamiltonian for the corresponding number of plane 
waves. A comparison of LCBB results with the conven- 
tional supercell approach was presented in Ref. Here 
we perform a convergence test, which consists in calculat- 
ing the energy levels with four different bulk basis sets of 
increasing size. Selected results are shown in Tab. | for 



3 



(GaAs)io/(AlAs)io and (GaAs)6/(AlAs)6 superlattices. 
As far as valence states are concerned, dependence of the 
energy levels on the basis set is below 10"'^ eV; for the 
lowest conduction band levels the dependence on the ba- 
sis set is generally below 0.05 eV, and falls below 10^^ eV 
when the 5*'' band of AlAs is included in the basis. The 
results of Tab .|| justify the use of basis 3), namely nt, = 1 
to 8 for GaAs and ni, = 5 for AlAs, as stated in Sec. II. 

In Figs. ^ and |^ we show the superlattice energy bands 
for n = 10: the electron energy levels are plotted along 
the highest symmetry lines in the tetragonal Brillouin 
zone (see Fig. |^). Since 2n k points in the fee Bril- 
louin zone are always folded onto the same k point in 
the smaller tetragonal Brillouin zone, the number of oc- 
cupied superlattice bands is 2n times the number of bulk 
bands for GaAs/ AlAs, n times the number of bulk bands 
for GaAs/vacuum superlattices. The dispersion along 
the r — Af direction is similar to the dispersion along 
the T — X direction in the bulk, while the other two di- 
rections r — .R and T — X have no counterpart in the 
band structures of Fig. ^. In GaAs/ AlAs, the dispersion 
along the growth direction T-Z is much smaller than in 
the other directions, as expected for superlattice mini- 
bands; in GaAs/vacuum the bands along T-Z are flat as 
tunneling through the vacuum has a negligible effect. 

The main differences in the superlattice band struc- 
tures compared to the bulk can be interpreted in terms 
of zone folding and quantum confinement effects; it is 
also interesting to compare the band structures of the 
two superlattices. Superlattice gaps are larger than bulk 
gaps: in particular GaAs/vacuum gaps are larger than 
GaAs/AlAs ones, as a result of a stronger confinement; 
moreover superlattice band gap widths increase as the su- 
perlattice period decreases. The lowering in the crystal 
symmetry is responsible for the removal of level degen- 
eracies: as an example in the GaAs/AlAs D2d superlat- 
tice the threefold degenerate valence states at F (spin- 
orbit is neglected) are split in a twofold-degenerate and 
a non-degenerate state, while in the GaAs/vacuum C2v 
superlattices the degeneracy is completely removed. 

In GaAs/vacuum bands we clearly see the appearance 
of states lying in the forbidden energy gaps. The lowest 
one lies in the gap from -10 to -6 eV, while two other 
ones lie in the optical gap from to about 2 eV. A fourth 
state can be recognized at -5 eV around the M point, 
while in other regions of the Brillouin zone it resonates 
with the energy bands. Indeed, four surface states or res- 
onances are expected from the presence of two dangling 
bonds at the two interfaces of each GaAs layers. We can 
identify the surface states by studying the behavior of 
the probability to find an electron along the growth 
direction z, averaged over the in-plane x, y coordinates. 
Taking as an example the conduction miniband states at 
F, where the potential profile is characterized by 0.5 eV 
deep wells in GaAs layers, we observe (see Fig ^ that 
an electron in a surface state (j = 39 in the exemplified 
case) has a high probability to be localized on the sur- 
face and a decaying probability to enter the GaAs layer; 



on the other hand an electron in a bulk state {j = 36 in 
the figure) has an oscillating probability to be found in 
the GaAs layers. Both states are obviously evanescent in 
vacuum. From the position of the uppermost occupied 
band (j — 40) it follows that GaAs/vacuum superlat- 
tices have a metallic behavior: this is an artifact of the 
calculations, as our aim is to simulate insulating multi- 
layers made of GaAs and a wide-gap oxide. We get rid 
of the problem by excluding surface states as initial or 
final states in interband transitions. This simulates the 
formation of interface states or defects in a GaAs/oxide 
superlattice, which would saturate the dangling bonds. 
The surface resonance cannot be easily eliminated, but 
it produces small efii'ects on the spectra, since it lies deep 
in the valence band. 



B. Optical spectra 

The calculated imaginary part of the dielectric function 
£2 (w) for both GaAs/AlAs and GaAs/vacuum systems 
are shown in Fig. |^ for different superlattice periods n. 
Here we average over the three orthogonal polarization 
directions to obtain a scalar dielectric function. The elec- 
tronic bands of GaAs and AlAs are similar and therefore 
the two optical spectra present the same features, namely 
El and E2 peaks (see Fig. ||). Superlattice spectra show 
also these structures, which can be compared directly to 
bulk spectra: in fact, even if in a superlattice selection 
rules allow more transitions as a consequence of zone- 
folding, the tetragonal Brillouin zone is smaller and, at 
last, all transitions which contribute to superlattice peaks 
have their equivalent counterparts in the bulk Brillouin 
zone. As a general remark, the Ei transition is found to 
blue shift and to split into two peaks for decreasing su- 
perlattice period; the confinement-induced shift is larger 
for the GaAs/vacuum system. On the other hand the E2 
transition is split for large-period GaAs/ AlAs superlat- 
tices, where the electronic states are confined in the two 
bulk layers leading to a superposition of the two bulk 
spectra; the two peaks merge into a single one for small 
period. A single E2 peak with a small blue shift is found 
for GaAs/vacuum superlattices. 

In Fig. H the peak energies are plotted as a function of 
the superlattice period n. First we comment on the be- 
havior of El : in zinc-blende crystals it comes from tran- 
sitions along the A line, in a region where bulk bands are 
almost parallel. When the system is confined in the [001] 
direction, it is not intuitive to describe the consequences 
of folding of (111) directions. The calculated spectra 
show that along the folded A line transitions subdivide in 
two main groups and lead to a splitting of the Ei peak in 
the absorption curves. The two peaks have different os- 
cillator strengths and, except for an intermediate period 
length, the lowest energy one becomes much stronger and 
covers the other one. Both peaks undergo confinement 
effects and are moved towards higher energies in com- 
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parison with their bulk position: confinement and the 
consequent shifts are stronger at smaUer well widths. A 
splitting of the Ei transition with a blue shift of bo 
peaks was indeed observed in GaAs/AlAs superlattice: 
In the present calculation this is attributed to a splitting 
of the bulk valence band at the L point and along the 
A line, as indicated by the band energies. The results 
of Figs. and ^ show also that Ei peak displacements 
are more relevant in GaAs/ vacuum superlattices, where 
quantum confinement effects are stronger due to the vac- 
uum barrier. 

The behavior of the E2 peak is substantially different: 
its main contribution comes from transitions in a region 
close to the special point k=^(|,i,i)j23At this point 
the alignment of both valence and conduction GaAs and 
AlAs bands is almost flat and the electronic wave func- 
tions are completely delocalized all over the heterostruc- 
ture. This explains why E2 peak in GaAs/ AlAs superlat- 
tices is at an intermediate energy between bulk GaAs and 
AlAs E2 peak positions and does not shift when the su- 
perlattice period n decreases. Our calculated peak 
tions are in good agreement with experimental dataj^ 
in particular a splitting of E2 is reported in Ref. ^ In 
GaAs/vacuum superlattices the situation changes: elec- 
trons near the special point are confined in GaAs lay- 
ers and the superlattice E2 peak has a weak blue shift 
at small superlattice periods n, going back to the bulk 
GaAs E2 energy when n grows. A single peak obviously 
arises in this case since there is no AlAs contribution. 

At last, we present in Fig. || some curves for the real 
part of dielectric function ei for GaAs/ AlAs heterostruc- 
tures. We observe that the average, or Penn gap (de- 
fined as the energy at which ei goes through zero) does 
not depend on the superlattice period. This proves that 
the center of gravity of valence and conduction bands 
is preserved, as suggested in Ref. |l^: this follows from 
compensating effects of a blue shift at the bottom of the 
band (positive curvature) and a red shift at the top of 
the band (negative curvature). 



C. Birefringence and absorption anisotropy 

In (GaAs)„/(AlAs)„ (001) superlattices the Td point 
group of the zinc-blende structure is replaced by the D2d 
symmetry group. Cubic crystals present isotropic optical 
properties, while in a superlattice the reduction in sym- 
metry leads to optical anisotropy in the real part of the 
dielectric constant (birefringence) and in the imaginary 
part (absorption anisotropy, or dicroism). The system is 
uniaxial, with the optical axis directed along the growth 
direction z, thus the dielectric tensor has the form: 

eij (uj) = en (uj) Sij , (3) 

where Cxx = ^yy — e_L and e^z — e||S The macroscopic 
dielectric tensor including local-field effects should be 
found by the procedure of inverting a matrix e(G, G') 



with off-diagonal jelepajents depending on the reciprocal 
lattice vectors G£2rEJ In bulk semiconductors the off- 
diagonal terms have a small effect on optical proper- 
ties, in superlattices however they are expected to give a 
sizeable-contribution and to depend on the superlattice 
period. E3 Nevertheless, at this stage we do not include 
local field effects in our calculation: we directly take into 
account only the effects of electronic confinement and 
band folding on optical transitions. Comparison with 



birefringence data reported in Ref. 16 should allow to 
determine if one of the two contributions to anisotropy 
is dominant. 

In Fig. |l^ we present the frequency dispersion of e|| 
and e_L, both for real and imaginary parts: the birefrin- 
gence is dispersionless up to energies close to the direct 
gap, while at higher energies it presents resonant contri- 
butions. We see as expected that folding and confinement 
can have a greater influence on the resonant part of the 
birefringence: indeed, transitions from valence subbands 
couple differently with xy or z polarized electric fields. 
Note that the interband absorption edge is higher in en- 
ergy for z-polarized light: this is in agreement with well 
known quantum well and superlattice physics, in which 
the lowest transition is a heavy, hole one and is forbid- 
den for light polarized along z.EJ Once again the effect is 
greater when confinement is stronger (small superlattice 
period). There is a dispersionless contribution to bire- 
fringence at low frequencies of the order of 10^'^ — 10^^ 
that cannot be distinguished in Fig. 10, As proposed in 
Ref. we decouple this low energy background birefrin- 
gence, describing An (uj) = {e±)^^^ — i^wY^^ terms of 
a fitting function 



An (uj) — Aniig — At 



gap 



In 1 - - 



(4) 



where Auhg is the background contribution we want to 
isolate, the second term refers to the resonant contribu- 
tion and ujg is the gap frequency. In Ref. |l^ the three 
parameters Arihg, Augap and ujg are extracted by fitting 
with expression (^) the experimental data. We also fit 
our calculated An (uj) curves by means of expression (jj). 
In Fig. ll| we display the fit parameters as a function of 
the well width: the graphs can be easily compared with 
the analogous experimental curves presented in Ref. |l^. 
The gap frequencies we extract by the fit agree both with 
directly calculated [gaps (barely visible in Fig.0) and with 
the measured ones.112 An^^p shows an increase of the reso- 
nance for small periods: the theoretical curve reproduces 
the trend of the experimental curve, although the calcu- 
lated values are smaller. The sudden decrease of the mea- 
sured background birefringence Antg below 40 A, on the 
other hand, is completely missing in our results. More- 
over the calculated magnitude of this term above 40 A 
is remarkably underestimated. The fact suggests that 
folding effects give only a minor contribution, while the 
origin of the behavior of An^g must be attributed to local 
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field effects, as already suggested in Ref. |16|. Similar con- 
siderations can be made for the GaAs/vacuum system. 

The magnitude of the background birefringence related 
to local field effects can be easily estimated in the case of 
long wavelength of incident radiation and not too small 
superlattice periods n. If boundaries are assumed abrupt 
and the constituent materials are supposed to conserve 
their bulk dielectric functions up to the interfaces, we can 
apply approximate expressions for e±_ and e|| in term s o f 
bulk constituent scalar dielectric functions, ei and e2t£l: 

e^(^)= , , {ei (uj) h + 62 {uj) h) , (5) 
n + h 

where h, I2 are the layer thicknesses of the two differ- 
ent materials. Using our calculated values of bulk static 
dielectric constants in expressions (|) and d) we find a 
rough estimate An ~ 0.03 for the local field contribution 
to birefringence. This value is much larger than the cal- 
culated values in Fig. |ll|, for all superlattice periods n, 
and it is of the same order of magnitude as the the exper- 
imental results at intermediate n. Only for small periods 
n, when relations (|^) and (||) do not hold, the local field 
correction should become small. Precise calculations of 
the optical properties including local-field effects are ob- 
viously required to clarify this point. 

IV. SUMMARY AND DISCUSSION 

The energy bands and optical constants of (001)- 
oriented (GaAs)„/(AlAs)„ and (GaAs)„/(vacuum)„ su- 
perlattices with n from 4 to 20 have been calculated by 
the LCBB method introduced in Ref. ^ This approach, 
in which the electronic states of the superlattice are ex- 
panded in the basis of bulk states calculated by empirical 
pseudopotentials, is found to be adequate and practical 
for superlattices with small to intermediate period; it is 
particularly useful for calculating how the optical spectra 
of the bulk materials are modified upon confinement. 

Quantum-confined induced shifts of the critical point 
energies are calculated for both kinds of superlattices 
and are found to be larger for the GaAs/vacuum system, 
where coupling between different GaAs layers is only due 
to quantum-mechanical tunneling and has a negligible ef- 
fect. For both GaAs/AlAs and GaAs/vacuum superlat- 
tices, the El peak in the absorption spectrum splits into 
two peaks with different blue shifts for decreasing super- 
lattice period. This result agrees with the observations of 
Ref. ^ on GaAs/AlAs superlattices, and is attributed to a 
symmetry splitting of the valence bands at the L point. 
The blue shifts are again larger for the GaAs/vacuum 
system. The E2 transition instead is found to be split 
for large-period GaAs/AlAs superlattices, where the elec- 
tronic states of the bulk are confined in each layer and 



the absorption spectrum is the superposition of the two 
bulk ones; the energy of the E2 peak depends weakly on 
the superlattice period. The average or Penn gap, de- 
fined as the first zero of the real part of the dielectric 
constant, does not depend on superlattice period, con- 
firming the expectation that a blue shift at the lower ab- 
sorption edges is compensated by red shifts in the upper 
parts of the absorption spectrum. The band contribution 
to linear birefringence of GaAs/AlAs superlattices is cal- 
culated and compared with recent experimental results 
of Ref. ^ The zero-frequency birefringence is found 
to be much smaller than the experimental results: the 
observed static birefringence-is attributed to local-field 
effects as already suggested.t£l The frequency-dependent 
part of the birefringence arising from band folding and 
quantum confinement increases with decreasing super- 
lattice period as found in the experiment, although the 
calculated values are smaller. 

The present work can be extended in several directions. 
Within the LCBB scheme relying on empirical pseudopo- 
tentials, a more precise calculation of the optical prop- 
erties requires the inclusion of spin-orbit interaction and 
of local-field effects: the latter are particularly important 
to account for birefringence of the superlattice. A com- 
plete description of the Ei peak requires the inclusion of 
excitonic effects, which account for half of the oscillator 
strength of the transition in the bulk. For superlattices 
with very short period, the expansion over bulk states be- 
comes unreliable (unless using a very large basis) and an 
ab-initio approach is required: it would be interesting to 
compare the results of ab-initio and LCBB calculations 
for superlattices with a period n from 4 to 8, say, where 
ab-initio methods are feasible and the LCBB approach is 
still meaningful. 



ACKNOWLEDGEMENTS 

The authors are grateful to M. Cardona, R. Del Sole, 
A. Stella and P. Tognini for several helpful suggestions 
and comments and to L. Reining and N. Vast for a critical 
reading of the manuscript. 



M.L. Clielikowsky and J.R. Cohen, Electronic Structure 

and Optical Properties of Semiconductors (Springer, Berlin, 

Heidelberg 1989). 
^ P.Y.Yu and M. Cardona, Fundamentals of Semiconductors 

(Springer & Verlag, Berlin, Heidelberg 1996). 
^ M. Garriga, M. Cardona, N.E. Christensen, P. Lauten- 

schlager, T. Isu, and K. Ploog, Phys. Rev. B 36, 3254 

(1987). 

* M. Alouani, S. Gopalan, M. Garriga, and N.E. Christensen, 
Phys. Rev. Lett. 61, 1643 (1988). 



6 



^ p. Tognini, L.C. Andreani, M. Geddo, and A. Stella, Phys. 

Rev. B 53, 6992 (1996). 
® C.W. Teng, J.F. Muth, R.M. Kolbas, K.M. Hassan, A.K. 

Sharma, A. Kvit, and J. Narayan, Appl. Phys. Lett. 76, 43 

(2000). 

J.-M. Gerard and B. Gayral, in Confined Photon Systems, 
edited by H. Benisty, J.-M. Gerard, R. Houdre, J. Rarity 
and C. Weisbuch (Springer, Berlin, 1999), p. 331; J.-M. 
Gerard, private communication. 

* A. Baldereschi, S. Baroni, and R. Resta, Phys. Rev. Lett. 
61, 734 (1988). 

^ The period of the superlattice is related to the number n 
of monolayers of each material by d = 2n5, where 5 is the 
distance between neighboring planes of the same atom. For 
(OOl)-grown superlattices, as considered here, d = na. 
L.J. Sham and M. Schliiter, Phys. Rev. Lett. 51, 1888 
(1983); Phys. Rev. B 32, 3883 (1985); J. P. Perdew and 
M. Levy, Phys. Rev. Lett. 51, 1884 (1983). 
" M. Hybertsen and S.G. Louie, Phys. Rev. Lett. 55, 1418 
(1985); Phys. Rev. B 34, 5390 (1986); R.W. Godby, M. 
Schliiter, and L.J. Sham, Phys. Rev. Lett. 56, 2415 (1986); 
Phys. Rev. B 35, 4170 (1987). 

S. Albrecht, L. Reining, R. Del Sole and G. Onida, Phys. 
Rev. Lett. 80, 4510 (1998); M. Rohlfing and S.G. Louie, 
Phys. Rev. Lett. 81, 2312 (1998). 

L.W. Wang, A. Franceschetti, and A. Zunger, Phys. Rev. 
Lett. 78, 2819 (1997). 
" L.W. Wang and A. Zunger, Phys. Rev. B 59, 15806 (1999). 
A. Fainstein, P. Etchegoin, P.V. Santos, M. Cardona, and 
K. Totemeyer, Phys. Rev. B 50, 11850 (1994). 
A.A. Sirenko, P. Etchegoin, A. Fainstein, K. Eberl, and M. 
Cardona, Phys. Rev. B 60, 8253 (1999). 
A. Fiore, V. Berger, E. Rosencher, P. Bravetti, and J. Na- 
gle. Nature 391, 463 (1988). 

This effect should be distinguished from circular birefrin- 
gence, or girotropy, which is also allowed by superlattice 
symmetry but is a fc-linear effect, i.e., an effect of spa- 
tial dispersion. See e.g. Landau-Lifshitz, Electrodynamics 
of Continuous Media (MIR, Moscow, 1969). 
V.M. Agranovich, and V.E. Kravtsov, Solid State Com- 
mun. 55, 85 (1985). 

M. Cardona, in Semiconductor Superlattices and Interfaces, 
Proceedings of the International School of Physics "E. 
Fermi", Course CXVII, edited by A. Stella and L. Miglio 
(North-HoUand, Amsterdam, 1993), p. 435. 
S.L. Adler, Phys. Rev. 126, 413 (1962). 

22 N. Wiser, Phys. Rev. 129, 62 (1963). 

23 W. Hanke, Advances in Physics 27, 287 (1978). 
2^ V.U. Nazarov, Phys. Rev. B 49, 17342 (1994). 

2® For a description of symmetry properties of superlattices 
see L.J. Sham and Y.-T. Lu, J. Lumin. 44, 207 (1989). 

2'' K.A. Mader and A. Zunger, Phys. Rev. B 50, 17393 (1994). 

2^ F. Bassani and G. Pastori Parravicini, Electronic States 
and Optical Transitions in Solids (Pergamon Press, Ox- 
ford, 1975). 

2* D.J. Chadi and M.L. Cohen, Phys. Rev. B 8, 5747 (1973). 

2^ S. Froyen, Phys. Rev. B 39, 3168 (1989). 

3° What we call E2 peak may in fact contain weak contribu- 
tions from the transition, which can be clearly resolved 
only in derivative spectra. See e.g. Refs. 



3^ E.E. Mendez, L.L. Chang, G. Landgren, R. Ludeke, L. 
Esaki, and F.H. PoUak, Phys. Rev. Lett. 46, 1230 (1981). 

■^2 In GaAs/vacuum superlattices, where the point group is 
C2v, the system is biaxial with principal axes along [110], 
[lIO], [001] and the dielectric constant for in-plane polar- 
ization has a slight additional anisotropy. The situation 
is analogous to heterostructures with no-common atom, 
where the in-plane anisotropy has been measured: see e.g. 
O. Krebs and P. Voisin, Phys. Rev. Lett. 77, 1829 (1996); 
O. Krebs, P. Voisin, D. Rondi, J.L. Centner, L. Goldstein, 
and J.C. Harmand, Physica E 2, 59 (1998). 

■^■^ The macroscopic treatment remains valid as long as the 
superlattice period remains much smaller than the wave- 
length of light, otherwise the superlattice should better be 
viewed as a one-dimensional photonic crystal, where local- 
field effects embodied in e(G, G') lead to the formation of 
a band gap for light propagation. 

3* The polarization selection rule follows immediately from 
k • p theory, although a precise calculation of absorption 
spectra close to the fundamental gap should of course in- 
clude the spin-orbit interaction. 




FIG. 1. Brillouin zone for simple tetragonal 

(GaAs),i/(AlAs)„ and (GaAs)n/(vacuum)„ (001) superlat- 
tices, included in bulk conventional cubic cell. The figure 
shows high symmetry points positions. 
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FIG. 2. Bulk band structure of (a) GaAs and (b) AlAs 
crystals along the high symmetry directions, obtained by em- 
pirical pseudopotentials in Ref. KfJ. The energy zero is always 
taken at the valence-band maximum of bulk GaAs and the 
valence band offset between GaAs and AlAs is 0.5 eV. 
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FIG. 4. Band structure of a (GaAs)io/(AlAs)io (001) su- 
perlattice along the high symmetry directions. The length of 
the T — Z line is multiplied by ten for clarity. The energy zero 
is taken at the bulk GaAs valence-band maximum. 
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FIG. 3. Calculated dielectric function (real and imaginary 
parts) of (a) GaAs and (b) AlAs. 



FIG. 5. Band structure of a (GaAs)io/(vacuum)io (001) 
superlattice along the high symmetry directions. The energy 
zero is taken at the bulk GaAs valence-band maximum. The 
uppermost occupied band is number 40 (dotted). 
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FIG. 6. Planar averaged probability along the growth di- 
rection z to find an electron at F in the 36th and the 39th 
band of a (GaAs)io/ (vacuum) lo (001) superlattice. 



FIG. 7. Imaginary part of the dielectric function for 
(GaAs)„/(AlAs)„ and (GaAs)„/(vacuum)„ superlattices, for 

different values of the period n. Different curves are offset 
for clarity. Ej splitting cannnot be easily seen in the figure: 
the peak positions have been determined by an enlargement 
of the spectral region of interest. 
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FIG. 8. El and E2 peak positions for GaAs/AlAs (closed 




FIG. 9. Real part of dielectric function for 

(GaAs)n/(AlAs)n (001) superlattices, for different values of 
the superlattice period n. 
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FIG. 10. (a) Components of the dielectric tensor (real 
and imaginary parts) for a (GaAs)4/(AlAs)4 (001) superlat- 
tice and (b) linear birefringence An = (ex)^^^ — (ey)^'^^ for 
(GaAs)„/(AlAs)„ (001) superlattices. 



FIG. 11. Parameters of equation (^) as a function of 
well width (half of a superlattice period): (a) energy gap, 
(b) Angap and Antg. The results can be compared with the 
analogous experimental graphs in Ref. |l6[ 



TABLE I. Comparison of energy levels (in eV) at the symmetry points for (GaAs)„/(AlAs)„ superlattices with period 
n = 6 and n = 10. We show highest valence levels and lowest conduction levels: the energy zero is taken at the valence band 
maximum. Four different choices are considered for the basis set: 1) five GaAs bands for nj from 1 to 5, no AlAs bands; 2) 
eight GaAs bands for Uh from 1 to 8, no AlAs bands; 3) eight GaAs bands for rih from 1 to 8, the 5th band of AlAs; 4) eight 
GaAs bands for nj, from 1 to 8, four AlAs bands for nj, from 5 to 8. 



basis 


F 


valence band maximum 
Z R M 


X 


F 


conduction band 
Z R 


minimum 

M 


X 


n=10 






















1) 


0.0000 


-0.0004 


-0.8992 


-2.1127 


-0.8992 


1.8318 


1.8303 


2.0069 


1.9048 


2.0063 


2) 


0.0000 


-0.0004 


-0.8990 


-2.1129 


-0.8990 


1.8173 


1.8190 


2.0057 


1.9014 


2.0049 


3) 


0.0000 


-0.0004 


-0.8990 


-2.1131 


-0.8990 


1.7949 


1.8008 


1.9911 


1.8588 


1.9881 


4) 


0.0000 


-0.0003 


-0.8985 


-2.1130 


-0.8985 


1.7884 


1.7957 


1.9904 


1.8590 


1.9874 


n=6 
1) 
2) 
3) 
4) 


0.0000 
0.0000 
0.0000 
0.0000 


-0.0118 
-0.0115 
-0.0115 
-0.0113 


-0.8488 
-0.8481 
-0.8481 
-0.8468 


-2.0968 
-2.0970 
-2.0973 
-2.0967 


-0.8489 
-0.8482 
-0.8482 
-0.8468 


1.9766 
1.9526 
1.9062 
1.9042 


1.9940 
1.9567 
1.9109 
1.9098 


2.1730 
2.1694 
2.1451 
2.1440 


1.9961 
1.9940 
1.9618 
1.9626 


2.1561 
2.1535 
2.1053 
2.1043 
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